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Abstract 

We describe the departure from equilibrium of matter distribu- 
tions representing sources for a class of Weyl metric. It is shown that, 
for extremely high gravitational fields, slight deviations from spherical 
symmetry may enhance the stability of the system weakening thereby 
its tendency to a catastrophic collapse. For critical values of sur- 
face gravitational potential, in contrast with the exactly spherically 
symmetric case, the speed of entering the collapse regime decreases 
substantially , at least for specific cases. 
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1 Introduction 



"L 'interviewer. - mais croyez-vous que vous devez etre d'accord 
avec ce que la plupart de gens autour de vous pensent?" 
"Ruth. - eh bien, c'est-a-dire que, quand je ne le suis pas, je me 
retrouve toujours a l'hopital..." 

R. D. Laing, A. Esterson, L'equilibre mental, la folie et la famille, 
Ed. Maspero (Pahs, 1971). 

As it is well known, since the early seminal work of Israel 0, the only 
static and asymptotically-flat vacuum space-time possessing a regular horizon 
is the Schwarzschild solution. For all the others Weyl exterior solutions 
0], the physical components of the Riemann tensor exhibit singularities at 
r = 2m. 

Since these components represent the true observables of the theory as it 
appears from the definition of tidal forces ||, it is intuitively clear that for 
high gravitational fields, the evolution of sources of Weyl space-time should 
drastically differ from the evolution of spherical sources It is important 
to keep in mind that the sharp difference in the behaviour of both types 
of sources (for very high gravitational fields) will exist independently on the 
magnitude of multipole moments (higher than monopole) of the Weyl source. 
This is so because, as the source approaches the horizon, any finite perturba- 
tion of the Schwarzschild space-time becomes fundamentally different from 
any Weyl solution, even when the latter is characterized by parameters whose 
values are arbitrarily close to those corresponding to the spherical symme- 
try. This point has been stressed long time ago ||, but usually it has been 
overlooked. 

It is the purpose of this work to study the evolution of axisymmetric 
sources for very high gravitational fields. This will allow us to put in evidence 
the role played by the non sphericity (however small) of the source, on the 
outcome of evolution. However, instead of following the evolution of the 
system long time after its departure from equilibrium, what would require 
the use of numerical procedures, we shall evaluate the source immediately 
after such departure. Here "immediately" means on a time scale smaller 
than hydrostatic time scale -see section 4 for more details. In doing so we 
shall avoid the introduction of numerical procedures. On the other hand, 
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however, we shall obtain only indications about the tendency of the object 
and not a complete description of its evolution. In spite of this limitation, 
this approach has proven to be useful in the study of spherically symmetric 
case (see and and references therein). 

As initial configurations, we shall consider two interior metrics. These 
ones were found some years ago by Stewart et al. || , following a prescription 
given by Hernandez |J allowing to obtain interior solutions of Weyl space 
time, from known spherically symmetric interior solutions. 

The configurations to be considered are sources of the so-called gamma 
metric (7-metric) ||10|, [II]. This metric, which is also known as Zipoy- 



Vorhees metric [Q, belongs to the family of Weyl's solutions, and is con- 
tinuously linked to the Schwarzschild space-time through one of its param- 
eters. The motivation for this choice is twofold. On one hand the exterior 
7-metric corresponds to a solution of the Laplace equation (in cylindrical co- 
ordinates) with the same singularity structure as the Schwarzschild solution 
(a line segment [K|). In this sense the 7-metric appears as the "natural" 
generalization of Schwarzschild space-time to the axisymmetric case. On the 
other hand, the two interior solutions considered have reasonable physical 
properties and generalize important and useful sources of the Schwarzschild 
space-time, namely the interior Schwarzschild solution (homogeneous den- 
sity) and the Adler solution JT3 . 



From all the comments above, it is not difficult to infer the astrophysi- 
cal relevance of the results here presented. Indeed, spherical symmetry is a 
common assumption in the study of compact self-gravitating objects (white 
dwarfs, neutron stars, black holes ) , furthermore in the specific case of non- 
rotating black holes, spherical symmetry should be " absolute" , according to 
Israel theorem. Therefore it is pertinent to ask, how do small deviations 
from this assumption, related to any kind of perturbation (e.g. fluctuations 
of the stellar matter, external perturbations, etc), affect the dynamics of the 
system?. The result obtained here, which deserves to be to emphasized, is 
that slight deviations from spherical symmetry seriously modify the depar- 
ture from equilibrium in the two examples presented. On the other hand, we 
are well aware of the fact that the 7 metric is not the only possible descrip- 
tion for the exterior of a compact objetc and, of course, the two equations 
of state considered here do not exhaust the list of possible candidates for 
the equation of state of the stellar matter. However, in view of the properties 
of the 7 metric and the two equations of state considered here, mentioned 
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above, it is fair to say that the case for the relevance of small deviations from 
spherical symmetry in the dynamics of compact objects has been established. 
This means that any conclusion on the structure and evolution of a compact 
object, derived on the assumption of spherical symmetry should be carefully 
checked against deviations from that assumption. 

This paper is organized as follows. In the next section we shall specify the 
space-time inside and outside the matter distribution, and the conventions 
used. In section 3 we give the energy momentum tensor components in 
terms of variables measured by a locally Minkowskian and comoving observer. 
The departure from equilibrium is analyzed in section 4, and results are 
discussed in the last section. We have included an appendix containing the 
components of the Einstein tensor, the Einstein equations in the limit of 
small non sphericity and the components of the conservation law. 



2 The Space-time 



2.1 The exterior space-time 

As has been mentioned above, our initial matter configuration is the source 
of an axially symmetric and static space-time (7- metric). In cylindrical co- 
ordinates, static axisymmetric solutions to Einstein equations are given by 
the Weyl metric 



with 
and 



ds ? = e 2X df 2 _ e -2A ^ (j p 2 + ^ + p 2 ^ 

A, pp + p -1 A, p + \ tZZ = 



fj, )g = 2p\ iP X y 



(1) 

(2) 
(3) 



Observe that (H) is just the Laplace equation for A (in the Euclidean space) 
The 7-metric is defined by [TO 

■Ri + R 2 -2m 



A = - In 



Ri + R 2 + 2m 



(Ri + R 2 + 2m) {R 1 + R 2 - 2m) 
AR X R 2 



(4) 
(5) 
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where 

Rl = p 2 + {z-mf R\ = p 2 + (z + mf . (6) 

It is worth noticing that A, as given by ([|), corresponds to the Newtonian 
potential of a line segment of mass density 7/2 and length 2m, symmetrically 
distributed along the z axis. The particular case 7 = 1, corresponds to the 
Schwarzschild metric. 



It will be useful to work in Erez- Rosen coordinates []P2"|, given by 

p 2 = (r 2 — 2mr) sin 2 9 z = (r — m) cos#, (7) 
which yields the line element as [ITU 



ds 2 = Fdt 2 - F- 1 {Gdr 2 + Hd9 2 + (r 2 - 2mr) sin 2 9d^ 2 ) , (8) 

where 



2m V 
r / 



2 o \ r 2 -! 

r — 2mr 



(9) 



G = — , (10) 

\ r 2 — 2mr + m 2 sin 9 I 



and 



h = - ^ I") 



(r 2 — 2mr) 

(r 2 — 2mr + m 2 sin 2 9^ 

Now, it is easy to check that 7 = 1 corresponds to the Schwarzschild metric. 

The total mass of the source is ]TU|, |TT[ M = 7m, and the quadrupole 
moment is given by 

Q = |m 3 (l - 7 2 ) • (12) 
So that 7 > 1 (7 < 1) corresponds to an oblate (prolate) spheroid. 

2.2 The interior space-time 

The metric within the matter distribution bounded by the surface 

r = a (13) 
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is given by 



9tt = f 1 

g rr = _ / 2(i- 7 ) A7 2 -2 S i- 7 2 



gee = -r*f^&^ 

3w = -r 2 / 2(1 ~ 7) sin 2 9 (14) 

where /, A, £ and $ are functions whose specific form depends on the model 
under consideration. 

The two cases to be considered here are the solutions reported in [§], 
namely 

1. The modified constant density Schwarzschild solution 



with 



and 



2 

A ( r ) = 1-^5 (16) 

4 

$M) = f + ^V^Wsin 2 ^ (18) 



V(r) = l + -(a-r), (19) 
a 



5 2 = — — , (20) 

where p ss denotes the energy density in the spherically symmetric limit 
(7 = 1) 

2. The modified Adler solution 

f{r) = A + Br 2 (21) 
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with 



and 



Cr 2 C 2 r 4 

E(r,0) = 1 + ^ + i7^sin 2 fl (23) 

(A + 3Br 2 ) 2/3 4{A + 3Br 2 f 3 

$(r,6) = (A + Br 2 ) + — V -±> S m 2 9 (24) 

V ' 4(A + 3Br 2 f 3 

(25) 

6 / 5m\ / m\ _1 , , , . 



.4 
B 



5t» 
2a 



2m 



1/2 



m 



2a 3 (l - ?f 



1/2 



2/3 



2m fl - 2*Y 
C = i ^ (27) 

« 3 (!-¥) 

Before closing this section, two remarks are in order: 

1. Since we are considering the source described in ( l4[) as an initial state, 
the time derivatives of functions /, A, E and $ will be in principle 
different from zero. 

2. Junction (Darmois) conditions are satisfied at the boundary r = a -see 
[Bfl for details. 



3 The energy momentum tensor 

In order to give physical meaning to the components of the energy momentum 
tensor in coordinates (t, r, 6, (p), we shall develop a procedure similar to that 
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used by Bondi [15j in his study of non static spherically symmetric sources. 
Thus, we introduce purely local Minkowski coordinates (r, x, y, z) defined 
by 



(It 



Pdt 



(28) 



dx 



yl-7 A -l+7 2 /2 S (l-7 2 )/2 rfr 



(29) 



(30) 



dz = rsin(0)/ "'dip. 



(31) 



Next, since we are assuming that our source does not dissipate energy, 
then the covariant components of the energy momentum tensor, as measured 
by a local Minkowskian and comoving with the fluid observer, will be 



T, 



[IV 



( p \ 

Pxx Pxy 

Pyx Pyy 

V o o o Pzz j 



(32) 



where p is the energy density and in general p xx ^ p yy ^ p zz and p xy = p yx . 
We may write (]32|) in the form 



T\w — (p "H Pzz) U^U V p Z zf]^iv~\~ \Pxx Pzz) k^k v -\- (p yy 
where r]^ u denotes the flat space-time metric and 

£/■„ = (!, 0, 0, 



% = ( 0, 1, 0, 



Pzz) I fJ'v ^^Pxyk^li/) ; 

(33) 



(34) 
(35) 



( o, 0, 1, 



(36) 



The components of the energy-momentum tensor T^ v in (£, r, <£>) coor- 
dinates are linked to ([33|) by 



(37) 



8 



where = dx u /dx fJ- is the local coordinate transformation matrix between 
Minkowskian coordinates and (t, r, 6, ip) coordinates and the Lorentz matri- 
ces are given by 

L\ = T L\ = U t = -Y Wi D i = Li = P i + {T - 1) ™ iU}i , (38) 
where 

w 2 =w 2 x + w 2 y r=-=L=, (39) 

and w x and w y denote, respectively, the velocity of a fluid element along the 
x and y (r and 9) directions, as measured by our local Minkowskian observer 
as defined by (p8|)-(|3T|). Observe that we are considering the case w z = 0, 



which means that the system preserves the reflection symmetry (besides the 
axial symmetry). 

The non vanishing components of AJ^ are 

A[ = P (40) 

A x = / l- 7A -l +7 V2 S (l-7 2 )/2 (41) 

A y = r/ 7(7-i) $ (i-7 2 )/2 (42) 

A; = rsin(tf)/ 1 ^. (43) 
Then, ( |3~T| ) readly gives 

T u = / 27 r 2 (p + p xx w 2 x + p yy w 2 y + 2p xy w x Wy) (44) 

(Tw x p + p xx w x A x + PyyWyA + p xy [w x A + w y A x ]) (45) 

(TWyP + P XX W X A + PyyWyAy + j9 X y [W y A + W X Ay] ) (46) 
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{T 2 w 2 x p + p xx A 2 x + PyyA 2 + 2 Pxy AA, 



T, 



re 



with 



So that 



r j(7-l) 2 ^-l+7 2 /2 ( | ) (l-7 2 )/2 S (l-7 2 )/2 x 

(T 2 w x w y p + A [p xx A x + PyyA y ] + p xy 



(V 2 w 2 y p + p xx A 2 + Pyy A 2 y + 2p xy AA l 
T w = r 2 / 2(1 " 7) sin 2 (6)p zz , 



A 2 + A T A. 



x^y 



A 



(r - 1) W X Wy 



A, 



A, = 1 + 
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(47) 

(48) 

(49) 
(50) 

(51) 
(52) 
(53) 



TfM/ {P + Pzz) UfJJ u Pzz9fJ,l>~\~ (Pxx Pzz) k^k v -\- (Pyy P l Z z) I fJ'V 4" 2p xy k(pl „) , 

(54) 

where U^, and 1^ are obtained after applying the boost velocity ( |3"8"D and 
the coordinate transformation (HO) -(EST) to (|3^-fl36l), 



(55) 



Up = T ( -«;J 1 -TA- 1 ^/2 S (l-7 2 )/2 ) _ Wj/r/ 7(7-l)$(l-7 2 )/2 ) Q ) 

fc M = ( -r^/T, / 1 ^A- 1 +^/2 S (l-7 2 )/2 Aa;; r/ 7(7-D$(l-7 2 )/2 A; Q 



/ M = ( -r^/T, / l-7 A -l+7 2 /2 S (l-7 2 )/2 A; r /7(7-l)$(l-7 2 )/2 A ^ Q ) . 



(56) 



(57) 
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4 Departure from equilibrium 



Let us consider a static axially symmetric source defined by (|TJ), which once 
submitted to perturbations, departs from equilibrium without dissipation. 
We shall evaluate the system after such departure, on a time scale such that 
w x and w y remain vanishingly small, whereas their time derivatives though 
small, will be different from zero. 

Thus, just after leaving the equilibrium, the following conditions hold 



Wy = W 



0. (i = r,6,(p) 
Wx,t, w y,t 7^ (small) 



(58) 
(59) 



From now on, unless otherwise stated, all equations are evaluated at the 
moment the system starts to deviate from equilibrium. 
Then from (^)-([50|), we obtain using (|58|) 



Tte — T tr — 
which implies, because of flT^D and (|75|) 

A, = f t = E t = $ t = 



(60) 
(61) 



where for simplicity we write for 0(u) (as we shall do hereafter). 

Obviously, spatial derivatives of the above quantities will be also vanish- 
ingly small on the time scale under consideration. 

Next, we shall evaluate the conservation law T£ = 0. 

The t-component yields 

P,t = (62) 
whereas the r and 9 component lead, after inspection of and Q85D, to 



f -^ f -^ A -l/2 A (e + eV2 )j: - { e + eV2 ){uJxAp + fe) + ^ ^ } 

2 

r 

r -lj-(2 £ + e 2 ) A -l/2 A ( £+£ 2 /2) < | ) (e+e 2 /2) E -( £ + £ 2 /2) 



L r(1 + e + a + _ t 



-2e - r 



Pxx Pxx, 



-^(-2e-e 2 ) + cot9 
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and 



f 



- + f -f(s) 
r f 



Pz 



+ 



i + ^( £ + e 2 ) + ^(-2 £ - £ - 



Pyy = 



j <P r 



x 



Pxy "I - Pxy,'i 



^l(-2e-e 2 )+cot9 



Pyy 



+ cot 6>p zz - + -^(-2e - e 2 )^* = 



(63) 



(64) 



where we have assumed 

7 = 1 + £ 

and e may be either positive or negative. 

Before going further into the analysis of the equations above it is quite 
instructive to consider the spherically symmetric situation. In this case we 
have 



E,0 — U)y }t — Pxy — P zz — Pyy — Pyy ft ~ S — 

Then, fl64j) becomes an identity and (|63"D reads 



(65) 



Vx,t{P +P* 



2 V.r , 

Pxx,r ~\Pxx Pyy) 2 \P P x 



>-5)/2 



where 



v r = 2 



m + A-Kr 3 p x . 
r(r — 2m) 



f = e 



u/2 



A = e' 5 = 1 - — 



2m 



with 



in ■ I \~r pdr 
lo 



(66) 

(67) 
(68) 

(69) 
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The physical meaning of fl66|) is quite transparent. It has the "Newtonian" 
form 

Force = mass x acceleration 

Indeed, the left-hand side consists of two factors. The time derivative of radial 
velocity and the inertial mass density. On the right-hand side, we have three 
possible sources of forces: The pressure gradient (negative), the possible 
anisotropic pressure contribution and the "gravitational" term (positive), 
multiplied by the relativistic correction factor e^ u ~ 5 ^ 2 . As the source becomes 
more and more compact the relativistic correction factor decreases as ~ r — 
2m, whereas the "gravitational" term grows as ~ 1/ (r— 2m). Since this latter 
term is positive and will prevail (as r — > 2m) over the two other force terms, 
we are lead unavoidably to a catastrophic collapse (u x ,t < 0), independently 
on the equation of state of the configuration. 

Let us now turn back to (|63|) and (|&iD, to infer what happens in the non- 
spherical case (i.e. e ^ 0), even though e « 1. Then neglecting higher order 
terms on e, and being careful with metric functions terms (some of which 
may tend to zero as the object becomes more and more compact) we obtain 
from flB3|) 

- /-[^ + »w)l=/A^(i)(, + fl .)(l + e) 

+ /A 1/2 A" £ S £ ^P xx>r - - [-2p XX + Pyy + P ZZ ] - E-^-(p xX ~ ^y) | 

- r- 1 / 1 " 2 ^ [cot 6p xy +p XVj e] (70) 
where we have used the fact that 

Vy, t ~ Pxy ~ Pzz - Pyy ~ 0(e) 

Since in this approximation the system (|63|), ( |64[ ) is not longer coupled (in 



u x j and uj y> t) we shall consider only equation (|70|). 

In order to extract more information from ( |7D| ) it is necessary to specify 
the source under consideration. We shall use the two configurations men- 
tioned in section 2. In both cases, / vanishes before a = 2m. Thus, in 
the Schwarzschild-like models we have /(0) = 0, if 2m/ a = 8/9. Since we 
know that any spherically symmetric static configuration with constant p ss 
and locally isotropic pressure should satisfy the constraint n = 2M/a < 8/9 
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we may assume that the system leaves the equilibrium for values of n close 
to 8/9. Then, an inspection of (|70|) shows that (if e > 0), for values of n 
approaching 8/9, the values of u x>t should become vanishingly small. Indeed, 
as n — > 8/9, / — > and p xx ~ 1/f ; /, r // ~ Pcca;> whereas A and E remain 
different from zero and bounded. Then, the largest term on the right of (ffO|) 
will be of order p xx (or 1/f). So, for any finite e (no matter how small) uj x t 
should be of order f 2e . 

We may now use Einstein equations to elucidate that for negative values of 
e the system may be unphysical. If n takes values close to the limit allowed by 
the model (8/9 for Schwarzschild-type model and 4/5 for Adler-type model), 
the critical values (/ — > 0, p xx — > 1/f ...) appear close to r = 0. Outside of 
this region the system is basically composed by an spherical incompressible 
fluid plus a perturbation in e. Thus, the physical or unphysical character of 
the model is determined by its behaviour close to these two limits. 

The energy density, in the limit r — > 0, for the Schwarzschild-type model 
and Adler-type model is given by expresions ( j8l|) and (f^) respectively. From 
these ones, it is easy to show that if e < 0, the energy density becomes 
negative as the system approaches to n — > 8/9 (Schwarzschild case) or n —>■ 
4/5 (Adler case) and positive for e > 0. Therefore, in both cases, positive 
energy conditions impose e > 0, and for very compact objects (close to the 
limit allowed by the model) the inertial mass term substantially increases 
suggesting the "stalling" of the collapse, close to, but before, the maximum 
allowed value of n. 



5 Conclusions 

We have seen so far that (as expected), for high gravitational fields, important 
differences appear between the spherical and the non-spherical collapse. This 
conclusion being true even for small non-sphericity. 

In the two models considered above, a factor multiplying the inertial mass 
term bring out those differences. Alternatively we may multiply both sides 
of equation ( ffO| ) by f 2e and say that the total "force term" decreases as f 2e . 
The final result being the same, namely that u X) t ~ f 2e - Although for the two 
examples considered here, we have considered e > 0, it is obvious that there 
exist models with e < 0, in which case the effective inertial mass density 
term may become very small, leading to highly unstable situations and to 
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a breakdown of the linear approximation. In fact, it is worth noticing that 
important differences between the two cases (e > and e < 0) appear also 
in the behaviour of the exterior 7-metric (|1U||, 

Thus, on the basis of presented results we may conclude that whatever 
the model and the sign of e would be, any source of Weyl metric would evolve 
quite differently from the corresponding spherical source, as critical values of 
n are considered. 

In particular, in the example examined, the inevitability of collapse in 
the spherical case, appears to be modified by a sharp increase in the effective 
inertial mass density term (or a sharp decrease in the "total force" term) as 
n approaches its maximum allowed value. This increase makes the system 
more stable, hindering its departure from equilibrium. 

We would like to conclude by stressing our main point: the inevitabil- 
ity of the catastrophic collapse for very compact objects, which appears 
in the spherically symmetric case, (for any regular matter configuration) 
is not present at least for the family of solutions considered here (e.g. is not 
"strictly" unavoidable). 
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A The Einstein tensor 

The calculation of non-vanishing components of the Einstein tensor 

Gfiu = Rpv ~ ~y$iivR (71) 

yields for the metric (0), using Maple V with GrTensor library and checking 
against the results given in [lJJ], 



G rr = -E 1 -^ 2 - 2 / 2-47 x 
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^(1 - 7) 2 + |f (1 - 7 2 ) + (^) 2 (7 4 - 1) + (^) 2 7(7 - 2)( 7 - l) 5 



2$ / 



(i- 7 2 )(2 7 2 -4 7 + i) 



-2 s l- 7 2 A 7 2 -2 y2-2 7 2 $7 2 -l 



l + |f(l-7 2 )cot£ 



4 + 1 



f r 1 9N 1 

+ l + 7 2 +- 



(72) 



Gw» = (1 - 7 2 



2E 7 + r + U"7 COt " 



(73) 



-^(7-l) 2 + |f(7 2 -l) 



+ (^(7 2 -2) + (l- 7 2 ) 
+ 1 f^(7 2 -2) + (l-7 2 ) 



St 


/,*" 


) 


— -\ 
2S 


/. 





^(l"7) 2 + ^(l-7 2 ) 



r V A 



2$ 



^(7 2 -l)(7-2) + |(l-7V^ 



(74) 



r 2 /27 2 -2 s7 2 -l A 2- 7 2 $ l- 7 2 x 



{ 

f \r f 

} 



^!i( 7 2_ 1 )_^:( 7 2_2) 

2E U ; 2A W ; 
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+r 2j2 7 ( 7 -2) $ l- 7 2 x 
{ 

^(7 2 -2) 



^(4-7 2 ) + ^(7 2 -l) + ^(4 7 -3) 



g(l+7 2 ) + ^(4 7 -3) 



+ f4l (l-7)(57- 



./ 



+ ^(2 -7 2 ) + t^(7 2 -1) 1) 



2A 



2E 



2/,t t , 
/ 



+ ^(l-7 2 )cot^ 



(l- 7 ^)cot^ + ^( 7 2 -l)-^(l-^)x 



^(l-7 2 ) + ^(7 2 -2)-^(l-7 2 )-^(l + 7 2 ) 



2E 



A 



2E 



{ 



2$ 



(7 ' - 11 (&< 72 + :) + k { '< 2 ' 2) - j* 27 ' 2 - 11 - ; 



^ : r / -i . \2/i \2 i ^,rr /-, 2\ , f,rr 2 

-2$2E (1 + 7) (1 ~ 7) + ^2$ (1 - 7) + T 7 



+ ^7 2 



2A 



+- 



^(2-7 2 ) + ^(7 2 -l) + ^(2 7 2 -l) 
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+r 2 sm 2 6f 2 -^ 



x 



{ 

24> 



+ 



./ 

fjt 
f 



^7(2 - 7 )( 7 - I) 2 + ( 7 2 -4 7 + 2) (^(2 - 7 2 ) - 



-^(l-7) 2 + ^(2 -7 2 ) + (7 2 -l) 



2E 2$ 



+ sin 2 ^ 2 ( 1 ^ 2 )$^ 2 - 1 x 



E 

~2E 



E^ 
2E 



(1-7)^ 



^(l + 7) 2 + ^(l-7 2 ) 



+^(l-7) 2 

{ 

^f(7 2 -2) 
2A U ; 



2$ 



(l-7 2 ) + ^(l"7) 2 



^(2 + 2 7 ) + ^(l-2 7 ) 



^(l-7) 2 + ^(l-7 2 ) 



y4 7 -2 A 2- 7 2 s7 2-l x 
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+ ^(1"7) 2 



^(l-7 2 ) + ^(l + 7) 2 
^(l-7 2 ) + ^7(l + 7) 



+ (1"7 4 ) ^ 



2$; 



+^(i 2 -i)- l f(i-i) 2 + 1 - 



^(7 2 -2) + |^(l-7 2 ) 



+- 



1 



£(3 t -1)(1- 7 )+3?£(y'-1)-± 



+r -2y2 7 (2- 7 ) $7 2-l x 
{ 

^( 7 2 -l) + ^ 



^(l-7) 2 (l + 7) 2 + §(l-7 4 ) 



2$ 2E 



cot + 1 



+ (l-7 2 ) 
} 

where a comma denote partial derivation. 



(78) 



B Energy density for 7 = 1 + £ 

Since we are interested in the effect of small deviations from spherical sym- 
metry, we shall assume that 

7 = l + e, (79) 

where e is a small constant that may be positive or negative. We restrict 
ourselves to first order in e, so we shall neglect terms of order 0(e 2 ) and 
higher. 

Then, the tt component of Einstein equations reads 



8vrp = y 2£ A £ S x 
{ 



19 



4 + £ 



1 

r 



i 2 



cot# 



S 



+ 1 



(80) 



B.l Schwarzschild-type solution 

Using ( p!5| - p0| ) into (|80| ) it is easy to show that, for r — > 0, the energy density 
can be expressed as 



P 



m 
2vra 3 



2c 



3 1 
- +£ 1 + 

2 I 3a /TT^_i 



(81) 



B.2 Adler-type solution 



In this case the energy density for r — > is given, by means of (21-^) and 
by expression 



3m / 1 — 5m/2a 

Ana 3 (1 - 2m/a) 1/3 Wl - 2m/a 



2e~2/3 



(l-m/a) 2/3 + 



1 - 5m/2a) 1/3 . 
(82) 



C Conservation equations 

The left hand side of conservation equations 



T" = 



are 



t& = r 2 x 
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£ (7 _ 1)(7 _ 2)+ ^_ 2)+ (£ + l) (1 _^) + ?£ x 

(P + P**^ + PyyW 2 y + 2p xy W x Wy) 

+P,t + Pxx,twl + 2p xx W x W x J + Pyy, t W 2 y 



+ 2PyyW y Wy ;t + 2p X y X W y + 2p X y (W X Wy >t + W X jW y) 

} 

+ j2 7 -l A l- 7 2 /2 S ( 7 2 -l)/2 r x 
{ 

(Tw x p + p xx w x A x + PyyWyA + p XJ/ [w x A + w y A x ]) 
+T, r w x p + Tw x , r p + Tw x p, r + p xx , r w x A x 

+Pyy, r WyA + ^ X y ^ [W X A + W <y A x ] 

+ WX,V [PXXA X + P X yA] + Wy, r \PyyA + J^Aj 

~\~A x ^ r [p xx w x 

+ PxyWy] + A, r [PyyWy + p xy W x ] 

} 



X 



X 



+r -lj2 7 -7 2 $ (7 2 -l)/2 r x 
{ 

iTWyP + P XX W X A + PyyWyAy + ^ X y [W y A + W X Ay]) 

+ T td Wyp + Wyfi (pT + PyyAy + P X y A) + W X fi {$> XX A + P X yAy) 
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+ A fi (PxxWx + PxyWy) + Ay td {PyyWy + p xy w x ) 



+ P, 9 TWy + PXX,0W X A + PyyfiWyky + ^ xyfi [iVyA + W X Ay] 

} 



+ 



^(7-D + i(l-7 2 ) 



f 



2$ 



X 



(^Wlp + PxX A 2 + P yyA 2 y + 2 Pxy AAy) 



+ f -fil-l)p z 



+ 



f " ' 2A V I ' 2S 
(r 2 ^p + Pxx Al + Pra A2 + 2p^AA x ) 



x 



_jl-2 7A -l+ 7 2 /2 S (l-7 2 )/2 r x 
{ 



/(t _ 1)(t _3) + |^- 2 ) + (| + |) (1 - 7 >) + ^ 

(r^P + P XX W X A X + PyyWyA + P X y [W X A + WyA x ]) 



+Tw x p it + T :t w x p + p xx>t w x A x + p yy ,tW y A + p xy>t [w x A + w y A x 



+ W Xtt (Tp + P XX A X + P xy A) + Wy <t (PyyA + P X yA x ) 



+A Xtt (p xy W y + p xx W x ) + A,t (PyyWy + p xy W x ) 



} 



f(T 2 -T + l) + |(l-f) + 



x 
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(T 2 w 2 x p + Pxx A 2 x + P yy h 2 + 2p xy hh x ) 



2IT r ty*p + 2T' z w x w Xtr p + T 2 w 2 x p >r + p XXy rK + P yy ,rA' z + 2p xy>r AA x 



-2A Xjr [p xx A x + p xy A] - 2A, r [p xy A x + p yy A] 



-ljl-7 2 A -l+7 2 /2 $ (7 2 -l)/2 S (l-7 2 )/2 



X 



{ 



^(l- 7 2 )+cot£ 
(r 2 w x w y p + p xx A x A + p yy A y A + p xy A 2 + A x A y ) 



+T 2 w x w y p fi + p xx ,eA x A + p yy fiA y A + p xyfi A 2 + A x A y 



+2TT fi w x w y p + V 2 w xfi w y p + V 2 w x w yfi p + A xfi [p xx A + p xy A y ] 



+A yfi [p yy A + p xy A x ] + h e [2p xy h + p xx h x + p^Ay] 

} 

-y7r 2 (p + PxxW* + p ra w 2 + 2p xy w x w y ) 



+ 



+ 



- + ^(1-7) 

1 , /,r 



Pz 



r + fT(7-l) + ^(l-T 2 ) 
(r 2 w 2 p + p ra A 2 + p ra A 2 + 2p xy hh v ) 



(84) 



{ 
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^ (7 _ 1) . + £( 7 .-2) + (^ + £)(l- 7 >) + E£ 

(YWyP + P XX W X A + PyyWyAy + X y [W y A + W X Ay]) 



X 



+ TWyP )t + PxX,tW X A + PyyjWyAy + P X y J [W y A + W X Ay 
+ T t fWyp + A >t [P XX W X + P X yWy] + AyJ [PyyWy + P X yW X 

+ w y,t [Pyy^y + PxyA + Tp] + W Xjt [p xx A + p xy A y ] 
} 



{ 



x 



(v 2 w x w y p + p xx A x A + Py y AyA + p xy A 2 + A x A y ) 



W X WyP J- -|- P XX ^A X A ~\- Pyy r AyA ~\~ P X y )T 



A 2 + A^ 



+2TT >r w x w y p + T 2 w x , r w y p + T 2 w x Wy yr p 
+A, r A x + p ra A y + 2p xy A] 

~^~Ay^ r \PyyA ~\~ P X yA X ] ~\~ A- X >r \j) XX A ~\~ J) ' X y Ajy] 
} 



Zj 1 

2^ (l - 7 2 ) + cot (T 2 w 2 y p + p ra A 2 + p ra A 2 + 2p xy AA y ) 



+Pzz COt 



24 



T 2 w 2 y p fi + Pxx,eA 2 + p yy , e A 2 y + 2p xyfi AA y 
2TT fi w 2 p + 2T 2 w y w yfi p 

-2k yfi [PyyAy + P xy A] ~ 2k fi [p xx A + Pxy^y] 

+ §! i 1 ~ 7 ') ( r ^ p + PxxA * + PyyA ' + 2p *y AA *) ^ 

= -p zz ^ (86) 
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